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ABSTRACT 

Simulation  is  a problem  solving  procedure  for 
defining  and  analyzing  a model  of  a system.  Computer-aided 
design  of  digital  logic  has  provided  the  design  engineer 
with  an  aid  to  reduce  the  tedious  and  time  consuming  task 
of  design  verification.  This  paper  describes  a simulation 
technique  for  the  simulation  of  digital  logic  circuits. 

This  paper  presents  a level  mode  logic  simulator 
that  has  improved  economy  in  execution  time  and  ease  of 
model  generation.  The  passage  of  time  is  simulated  in  a 
precise  fashion  and  element  models  are  executed  only  when 
activity  occurs  in  the  circuit.  A behavior  model  descrip- 
tion is  accomplished  on  an  element  level  rather  than  a gate 
level.  The  use  of  three-valued  logic  and  the  use  of  precise 
timing  delays  for  both  rising  and  falling  signal  levels 
present  a very  accurate  and  informative  circuit  output 
timing  diagram.  Th’s  is  demonstrated  by  simulation  of  an 
even-odd  detection  circuit  and  an  up-down  gray  code  counter. 


CHAPTER  1 


INTRODUCTION 

Simulation  is  a problem  solving  procedure  for  de- 
fining and  analyzing  a model  of  a system;  in  particular,  a 
digital  computer  simulation  is  the  establishment  of  a 
mathematical  or  logical  model  of  a system  and  the  manipula- 
tion of  that  model  on  a digital  computer.  Designers  of 
digital  logic  circuitry  have  long  recognized  the  need  for 
computer-aided  design  and  have  placed  much  emphasis  on 
digital  computer  simulation  of  digital  logic  circuits. 

Menial  and  time  consuming  tasks  such  as  logic  verification 
are  now  routinely  carried  out  by  computer  programs.  Several 
types  of  logic  simulators  ranging  from  simple  one-of-a-kind 
circuit  simulations  to  very  complex  programs  are  found 
throughout  government,  industry,  and  education  (Hill  and 
Peterson,  1974).  These  simulations  permit  inferences  to  be 
drawn  about  systems;  without  building  them,  if  they  are  only 
proposed  systems;  without  disturbing  them,  if  they  are 
operating  systems;  and  without  destroying  them,  if  the 
object  of  an  experiment  is  to  determine  their  limits  of 
stress . 

Without  the  assistance  of  computer  simulation, 
verification  of  the  design  of  any  large  digital  circuit  is 
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a manually  tedious  and  time  consuming  task.  The  design 
engineer  requires  no  creative  skills  to  accomplish  this 
routine  task.  Computer  simulation  relieves  the  designer  of 
this  task,  therefore  reducing  the  time  and  expense  involved 
in  testing  designs  before  fabrication  (Williams,  1975a, 
1975b).  Once  the  basic  circuit  operation  has  been  verified, 
the  designer  can  optimize  his  design  and  quickly  re-verify 
the  circuit  operation.  The  design  engineer  thus  has  a very 
powerful  tool  in  a simulator  to  assist  him  in  obtaining  a 
more  nearly  optimal  design  in  less  time. 

The  purpose  of  this  study  was  to  develop  a computer 
program  that  would  perform  a simulation  of  large  digital 
logic  circuits.  The  constraints  placed  on  the  program  were: 

1.  The  input  format  for  element  descriptions  was  to  be 
compatible  with  existing  programs. 

2.  The  simulator  would  be  a level  mode  simulator  where 
level  mode  means  that  a delay  time  is  associated 
with  each  logic  element.  It  should  be  noted  that 
this  constraint  assures  that  the  program  will  not  be 
restricted  to  a particular  family  of  logic  elements 
nor  a particular  type  of  logic  (synchronous, 
asynchronous,  combinational,  sequential). 

3.  Equivalent  models  for  devices  will  not  be  used.  For 
example,  AND  and  OR  logic  will  not  be  used  to 
describe  systems  of  NAND  or  NOR  logic,  or  flip-flops 


will  not  have  to  be  expanded  to  their  gate  level 
equivalence . 

4.  The  output  will  be  of  the  form  of  a timing  diagram. 
The  program,  called  Event  Simulator  (EVESIM),  is  primarily 
intended  for  use  by  the  students  and  faculty  of  the 
Department  of  Electrical  Engineering,  The  University  of 
Arizona. 

The  study  began  with  a thorough  review  of  available 
literature  on  digital  logic  simulation.  Chapter  2 briefly 
summarizes  previous  work  and  techniques  while  presenting  the 
basic  concepts  of  EVESIM. 

Structured  programming  techniques  have  been  applied 
and  basic  functions  performed  by  the  program  have  been 
implemented  with  separate  routines.  Using  these  separate 
routines  is  supported  in  three  advantages : 

1.  Debugging  of  the  program  is  simplified  since  all 
subroutines  could  be  checked  out  separately  before 
being  added  to  the  main  program. 

2.  With  minimal  effort,  someone  unfamiliar  with  the 
program  can  gain  familiarity. 

3.  Maintenance  and  improvements  to  the  programs  can  be 
done  easily  by  changing  or  adding  individual 
routines . 

The  use  of  separate  routines  for  separate  functions  is 
readily  recognizable  in  the  detailed  description  of  EVESIM 
as  presented  in  Chapter  4. 
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Since  it  is  beyond  the  scope  of  this  study  to 
develop  element  models  for  all  the  available  logic  compo- 
nents, Chapters  5,  6,  and  7 are  dedicated  to  providing 
sufficient  information  to  enable  a user  to  //rite  and 
implement  additional  element  models.  EVESIM,  as  it  exists, 
forms  a basic  system  that  is  designed  to  allow  easy  expan- 
sion. Future  work  to  enhance  the  capabilities  of  EVESIM 
might  include:  optional  data  inputs/outputs  and  a random 
delay  time  generator.  Further  recommendations  are  included 
in  Chapter  8. 


CHAPTER  2 


DIGITAL  SYSTEMS  SIMULATION 

The  usefulness  and  purpose  of  digital  logic  simula- 
tors has  been  well  established  in  Chapter  1.  It  is  now 
appropriate  to  review  the  many  approaches  to  logic  simula- 
tion. Hill  and  Peterson  (1974)  discuss  the  relation  of  cost 
to  useful  information  output  for  several  approaches  to  logic 
simulation.  The  cost  of  operation  of  simulators  for  large 
circuits  is  great,  thus  providing  the  incentive  to  develop 
efficiency  in  simulation  techniques. 

History  of  Simulation 

Logic  simulators  continue,  as  in  the  past,  to  be 
developed  based  on  how  much  information  is  required  about  a 
logic  circuit.  They  are  certainly  influenced  by  the  size  of 
the  circuit  and  type  of  circuit.  Some  simulators  are  part 
of  sophisticated  computer-aided  logic  design  packages  as 
developed  by  Hayes  (1969),  and  others  are  logic  verifiers  as 
written  by  Williaims  (1975a,  1975b).  Simulators,  however, 
normally  appear  at  least  once  in  the  automated  design 
process  prior  to  the  design  being  committed  to  hardware. 

In  the  case  of  Chappell,  Elmendorf , and  Schmidt  (1974)  , 
several  different  simulators  are  used  for  different  purposes 
in  the  design  process.  For  example,  using  a level  mode 


simulator  will  not  necessarily  insure  that  a design  is  valid 
over  the  allowed  range  of  element  delay  parameters.  The 
level  mode  simulator  can,  however,  be  used  to  predict  worst 
case  timing  conditions  and  to  provide  some  verification  to 
clock  mode  design  that  strict  adherence  to  the  clock  mode 
design  rules  has  been  accomplished.  The  use  of  a clock  mode 
simulator  as  presented  by  Williams  (1975a,  1975b)  for 
verification  of  the  combinational  logic  would  certainly 
enhance  a level  mode  simulator,  thus  providing  the  incentive 
for  developing  different  simulators  for  different  purposes 
in  the  design  process. 

Simulation  programs  differ  greatly  in  nature  and 
complexity  as  a result  of  the  variety  of  applications. 
Simulators  such  as  electronic  network  analysis  as  presented 
by  Anwaruddin  (1969)  are  very  complex  and  costly  to  operate, 
however  they  provide  the  most  accurate  simulation.  Simu- 
lators also  differ  in  regard  to  the  type  of  network  specifi- 
cations accepted  as  input.  For  example,  Stockwell  (1962) 
uses  Boolean  equations  for  circuit  descriptions  while  others 
use  the  actual  circuit  elements  and  their  interconnections, 
as  does  Williams  (1975a,  1975b) . Larson  and  Mano  (1965) 
have  presented  a simulation  technique  which  uses  appropriate 
digital  network  models  and  FORTRAN  IV  to  economically 
simulate  logic  networks  where  the  high  cost  of  a general- 
purpose  simulator  is  not  justified.  Stang  (1968)  and 
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Weingarten  (1965)  have  taken  similar  approaches  by  using  the 

logical  and  arithmetic  features  of  FORTRAN. 

Probably  the  area  of  greatest  interest  in  simulators 

is  level  mode  simulators.  Hill  and  Peterson  (1974)  discuss 

these  in  great  detail.  Again,  the  amount  of  information 

required  versus  cost  determine  the  approach  taken  in  a level 

mode  simulation.  In  order  to  make  detailed  analysis  of 

races  and  hazards  in  circuits,  an  analog  output  level  mode 

simulator  may  be  required.  This  type  of  simulation  is 

accomplished  by  usng  quantized  voltage  levels  between  the 

logic  values  of  zero  and  one.  Neglecting  rise  time  and 

only  delay  results  in  a 2-level  level  mode  simulation.  Such 

a simulator  using  an  up-down  counter  to  model  delay 

presently  exists  and  is  accessible  on  the  DEC  10  computer 

system  at  The  University  of  Arizona.  The  most  basic  ! 

I 

approach  to  simulation  is  a clock  mode  simulation.  One  j 

iteration  through  the  simulation  is  equivalent  to  one  clock 

pulse.  I 

I 

Great  effort  has  been  placed  in  development  of  I 

advanced  techniques  for  improving  the  calculation  and  effi-  j 

ciency  of  the  basic  simulation  models  to  be  discussed.  I 

: 

I , 

Event-directed  techniques  (Breuer,  1970)  and  exclusive  ;j 

simulation  of  activity  techniques  (Ulrich,  1969)  are  in-  j 

I 

tended  to  reduce  the  amount  of  calculations  required  in  one  j 


iteration  of  a simulation. 
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Modeling  Digital  Logic 

As  stated  previously,  simulation  is  a problem 
solving  procedure  for  defining  and  analyzing  a model  of  a 
system.  The  value  of  simulations,  of  course,  is  that  large 
amounts  of  data  about  a system  can  be  collected  in  a short 
period  of  time  under  controlled  conditions.  How  accurately 
the  behavior  of  the  real  system  is  depicted  by  the  data  is 
reflective  of  the  accuracy  of  the  simulation  model  and  the 
model  characteristics  (Szygenda  and  Thompson,  1975) . The 
simulator  designer  must  first  decide  whether  the  system  is 
to  be  used  for  logic  verification  or  design  verification. 

He  needs  to  ascertain  whether  to  use  two,  three,  or  multiple 
gate  output  values.  There  is  also  the  question  of  whether 
to  use  zero,  unit,  assignable,  or  precise  delays;  and 
closely  related,  whether  to  use  differing  signal  rise  and 
fall  times.  Certainly,  the  type  of  internal  device  model 
must  be  considered,  whether  it  be  gate  or  element. 

Logic  verification  would  be  similar  to  the  operation 
of  a combinational  logic  circuit  in  that  it  is  purely  a 
computation  of  all  possible  outputs  for  all  possible  input 
combinations.  It  is  also  equivalent  to  a clock  mode  simula- 
tion. Since  the  implementation  of  such  a simulator  is  very 
simple  and  if  no  other  information  is  required,  then  the 
simulator  designer  should  consider  only  logic  verification. 
If  more  information  is  desired,  such  as  timing  considera- 
tions, then  a more  complex  approach  is  appropriate,  noting 


that  lc)<|ic:  verification  can  be  a special  case  of  a more 
complex  simulator. 
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It  is  obvious  that  a logic  simulator  requires  at 
least  two  logic  levels,  0 and  1.  Using  only  two  values, 
however,  has  two  major  drawbacks.  The  first  disadvantage  is 
that  establishing  initial  values  on  all  elements  is  very 
difficult,  except  for  clock  mode  simulation.  For  large 
circuits,  establishing  initial  element  output  values  then 
becomes  costly.  The  second  disadvantage  is  that  unknown 
conditions  due  to  spikes  or  other  sources  can  not  be  repre- 
sented. The  third  value  normally  added  is  X or  unknown 
which  allows  simulators  to  overcome  the  disadvantages  of 
two-valued  logic.  Chappell  et  al.  (1974)  define  a non- 
propagating unknown  condition  that  is  used  during  the 
process  of  establishing  initial  values.  This  non- 
propagating feature  prevents  an  unknown  condition  from  over 
writing  an  already  established  known  value.  They  use  a ^ 

fourth  value  called  a propagating  unknown  for  the  remainder 
of  their  simulation.  It  is  feasible  to  expand  further;  for  i 

example  0,  1,  2 (non-propagating  unknown),  3 (propagating 
unknown) , 4 (signal  rising) , and  5 (signal  falling) . 

Certainly,  the  addition  of  4 and  5 would  be  useful  in 

simulating  rise  and  fall  times  if  that  information  was  j 

required . 

It  is  important  that  the  difference  between  gates 
and  elements  be  established.  A gate  is  defined  as  a device 

■ 
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that  has  multiple  inputs  and  a single  output.  An  element  is 
defined  as  a device  that  has  multiple  inputs  and  multiple 
outputs;  in  addition,  the  ordering  of  the  inputs  and  outputs 
is  significant.  A good  example  of  an  element  is  a flip- 
flop;  certainly  the  inputs  and  outputs  are  ordered.  These 
differences  result  in  major  differences  of  implementation. 

The  complexity  of  implementing  element  models  is  the  reason 
that  most  simulators  use  only  gate  models.  There  are  some 
serious  drawbacks,  however,  in  using  only  gate  models.  The 
first  of  these  is  that  a large  number  of  gates  are  required 
to  implement  some  elements,  for  example,  a JK  flip-flop  may 
require  eight  to  ten  gates  to  implement.  Another  major 

i 

drawback  is  that  the  interaction  of  gates  in  implementing  an  ! 

element  may  not  present  the  same  operating  characteristics  I 

I 

as  the  element  itself,  thus  presenting  a poor  model  of  the  i 

element.  Therefore,  careful  consideration  must  be  given  to  |j 

the  type  of  internal  device  models. 

To  achieve  greater  accuracy  in  the  model,  an 
assignable  or  precise  timing  capability  is  required.  Zero 

i 

and  unit  delay  (clock  mode)  simulators  are  only  useful  for  ! 

i 

logic  verification.  Timing  is  inherent  in  design  verifies-  j 

' i 

tion.  Assignable  delay  simulators  allow  assignment  of  a j 

I 

single  average  delay  to  each  element  or  gate  type.  For 
excimple,  one  AND  gate  in  a circuit  may  be  assigned  5 units 

of  delay  and  another  AND  gate  in  the  same  circuit  may  be  I 

i 

assigned  4 units  of  delay.  The  up-down  timing  technique  1 
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mentioned  earlier  is  a good  example.  This  model  is  made 
somewhat  more  precise  if  delay  assignments  are  made  on  an 
individual  element  based  on  such  parameters  as  fanout.  This 
technique  also  allows  circuits  consisting  of  different 
families  of  logic  to  be  accurately  simulated.  This  can  be 
further  expanded  to  assignment  of  different  delays  dependent 
upon  whether  the  signal  is  rising  from  logical  0 to  logical 
1,  or  falling  from  1 to  0 . Though  this  technique  is  not  a 
true  rise  and  fall  time  simulator,  as  needed  in  electronic 
analysis,  it  does  present  a more  accurate  picture  of  the 
real  world  situation. 

Implementation  of  Simulators 
Having  looked  at  modeling  techniques  in  great 
detail,  it  is  now  appropriate  to  answer  some  questions  about 
implementation  prior  to  discussing  EVESIM.  At  a minimum, 
areas  of  interest  in  implementation  should  include  the  ease 
of  use,  independence  of  data  processing  equipment,  selective 
tracing  and  time  flow  mechanisms,  and  element  descriptions. 

Engineers  who  make  extensive  use  of  computer  pro- 
grams, invariably  avoid  programs  that  are  difficult  to  use. 
Good  programs  minimize  user  required  tasks.  This  requires 
complex  pre-  and  post-processing,  but  must  be  done  in  order 
to  gain  wide  usage  of  the  progreim. 

Machine  independence  of  the  program  is  very  impor- 


tant, particularly  if  the  program  is  large  and  was  costly  to 
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develop.  A proqram  is  not  very  useful  if  the  author  of  the 
proqram  is  the  only  one  that  can  use  it  because  his  machine 
is  unique.  Though  low-level  languages  normally  provide  a 
savings  in  run  time,  they  are  not  very  transportable. 

Higher  level  languages  such  as  ANSI  standard  FORTRAN  have 
proven  to  be  much  more  transportable  and  therefore,  a better 
language  to  use. 

Gate  and  element  have  previously  been  defined;  how- 
ever, as  far  as  implementation  is  concerned,  a gate  is  a 
special  case  of  an  element.  Defining  the  operation  of  an 
element  in  terms  of  a subroutine  or  in-line  code  is  rela- 
tively straightforward  assuming  the  programmer  has  a through 
understanding  of  the  element's  operation.  Defining  an 
element  in  a subroutine  has  obvious  advantages  such  as  ease 
of  debugging,  ease  of  being  replaced  or  overwritten  by 
future  users,  and  ease  of  calling  for  use  throughout  the 
program.  A detailed  discussion  of  how  to  write  element  sub- 
routines is  contained  in  Chapter  5. 

Breuer  (1970)  has  suggested  a partitioning  of  a 
circuit  into  subsets  of  elements  U.(i)  such  that  if  after 
simulating  element  i,  it  is  found  that  i changes  state,  then 
simulate  the  elements  of  U^d)  where  the  elements  of  (i) 
are  those  elements  that  have  i as  an  input.  This  is 
selective  tracing;  in  other  words,  if  an  element  output  does 
not  change  when  new  inputs  have  been  evaluated,  then  the 
element  output  is  not  followed.  Since  partitioning 


alqorithms  are  very  complex,  Ulrich  (1969)  has  suggested 
that  this  partitioning  be  limited  to  the  fanout  of  a single 
element.  Breuer's  (1970)  technique  has  been  named  event- 
directed  simulation,  where  event  means  that  an  element 
output  has  changed.  The  simulator  jumps  from  event  to  event 
as  opposed  to  being  time  driven.  Ulrich's  (1969)  technique 
can  be  labeled  time-driven  event  simulator.  The  time-driven 
simulator  operates  by  incrementing  a fixed  time  unit  and 
checking  for  activity  at  that  time  where  activity  is  defined 
as  elements  scheduled  to  be  evaluated.  It  should  be  obvious 
that  these  techniques  are  required  for  large  circuits  since 
any  system  relying  on  evaluating  each  element  on  each 
increment  of  time  quickly  becomes  too  time  consuming  and  in 
turn  costly. 


Design  Concepts  of  EVESIM 

Having  established  a background  of  basic  simulation 
concepts,  it  is  now  time  to  begin  a more  detailed  analysis 
of  EVESIM.  Chapter  1 well  established  that  EVESIM  should  be 
a design  verification  simulator  and  that  it  should  use  an 
element  approach  to  modeling;  therefore,  these  two  subjects 
warrant  no  further  discussion  at  this  time.  EVESIM  uses 
three- valued  logic,  it  applies  precise  delay  techniques  that 
reflect  different  rise  and  fall  times,  and  it  uses  a hybrid 
of  event-directed  and  time-driven  event-directed  control 


techniques . 
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EVESIM  uses  three-valued  logic  and  in  particular,  i 

i 

the  unknown  condition  X,  for  two  primary  purposes.  The  J 

first  use  is  similar  to  that  of  the  non-propagating  unknown  j 

used  by  Chappell  et  al.  (1974).  During  the  process  of  - 

establishing  initial  values  for  the  elements  of  the  circuit 
to  be  simulated,  all  element  outputs  are  first  set  to  X, 

unknown,  except  for  memory  type  elements  which  are  user  ' 

specified.  Memory  element  outputs  are  not  allowed  to  change 
during  the  process  of  updating  the  other  elements.  In  a 
circular  process,  new  outputs  for  the  remaining  elements  are 
calculated  until  no  more  unknown  conditions  exist  in  the 
circuit.  A default  mode  and  an  override  capability  are 
built  into  the  system  in  case  all  unknown  conditions  cannot 
be  cleared.  (This  will  be  discussed  further  in  Chapter  4.) 

The  second  use  of  the  unknown  is  in  a propagating  mode  to 
reflect  illegal  conditions  such  as  timing  problems  generated 
in  the  circuit  being  simulated. 

Though  the  use  of  precise  delays  for  modeling 
normally  means  more  complex  programs,  precise  delays  were 
used  in  EVESIM.  Based  on  the  use  of  the  hybrid  driving 

technique,  to  Vje  discussed,  this  complexity  has  been  offset.  I 

EVESIM  allows  individual  element  assignment  of  differing 
rising  signal  delay  time  and  falling  signal  delay  time, 
thus  presenting  a more  accurate  model. 

The  heart  of  EVESIM  is  the  driving  mechanism,  as  it 
is  with  :-inv  simulation  program.  Ulrich  (1969)  has  reported 

I' 

J 


that  only  one  per  cent  of  all  elements  in  a typical  simu- 
lated digital  network  are  simultaneously  active.  Assiuning 
this  to  be  true,  then  any  driving  technique  that  takes 
advantage  of  this  is  an  improvement  over  conventional  tech- 
niques of  simulating  each  element  in  every  increment  of 
time.  This  suggests  a discrete  event  type  of  simulation. 

In  other  words,  it  is  desirable  to  simulate  an  element  only 
when  it  is  active.  Define  a set  C,  the  elements  of  C being 
the  elements,  including  external  inputs,  of  a circuit  C. 
Then,  at  any  instant  of  time  a subset  of  C,  consisting  of 
currently  active  elements,  is  the  only  source  of  future  net- 
work activity,  signals  originating  in  A^^  are  transmitted  to 
associated  destination  elements,  forming  a destination  set 
also  a subset  of  C.  Not  all  elements  of  respond  to 
the  stimulation,  thus  a new  subset,  A2 # of  active  elements 
is  formed  where  A2  is  a subset  of  EVESIM  takes  full 

advantage  of  this  concept. 


CHAPTER  3 


DRIVING  MECHANISM  OP  EVESIM 

Until  now,  simulator  driving  mechanisms  have  been 
discussed  in  a general  sense.  This  chapter  carries  the 
discussion  of  driving  mechanisms  to  that  used  in  EVESIM. 

The  information  required  for  the  modeling  of  elements  in 
the  driving  mechanism,  is  also  identified. 

The  driving  mechanism  used  by  EVESIM  is  best 
identified  as  being  a time-driven  event-directed  mechanism 
where  an  event  is  defined  as  an  element  output  taking  on  a 
new  value.  The  technique  is  a hybrid  of  a time-driven  and 
and  event-driven  mechanism.  Basically,  a time-driven 
circuit  model  is  incremented  by  a fixed  time  unit  and  a new 
set  of  element  output  values  is  calculated  per  time  incre- 
mented. An  event-directed  approach  is  driven  by  jumping 
from  event  to  event  irrespective  of  time.  It  should  be 
obvious  that  considerable  calculations  are  required  for  the 
time  driven  approach,  and  that  a considerable  amount  of  data 
must  be  maintained  for  each  event  in  the  event-directed 
approach.  The  technique  used  by  EVESIM  is  a combination  of 
these  two  approaches. 

EVESIM  is  time  driven  in  the  sense  that  events  can 
occur  only  at  points  in  time  that  are  an  integral  multiple 
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of  the  fixed  increment  of  simulation  time.  In  other  words, 
simulation  time  is  merely  represented  by  a counter,  TIME^^^ 

= TIME^+1,  beqinninq  at  time  0,  and  events  can  occur  only 
at  these  integer  values  of  TIME.  Element  models  are  used 
to  determine  new  output  values  as  element  inputs  change. 

The  model  must  also  determine  when  the  element  output  is  to 
change  in  terms  of  simulation  time.  This  timing  information 
is  calculated  by  the  model  from  user  input  element  delay 
times.  This  time  of  change  or  time  of  occurrence  for  an 
event  is  used  to  properly  insert  the  event  into  an  event 
list.  This  event  list  is  the  only  source  of  future  activity 
in  the  circuit  being  simulated.  In  other  terms,  these 
events  or  known  element  output  changes  are  the  only 
occurrences  in  the  circuit  that  can  cause  future  changes  of 

other  elements.  The  event  is  executed  by  following  the 

( 

element  output  to  all  fanout  elements  connected  and  checking 
the  fanout  elements  for  future  activity.  Therefore,  simu- 
lating only  the  effects  of  these  changes  has  reduced  the 
simulation  to  an  event  directed  simulation.  Allowing  these 
events  to  occur  only  at  fixed  increments  of  time  results  in 
a time  driven  event  directed  simulator. 

To  accomplish  the  tasks  of  driving  the  simulation, 
EVESIM  uses  two  arrays.  The  first  is  ICL(N,2)  where  N is 
equal  to  the  maximum  delay  time  assigned  to  an  element  of 
the  circuit  being  simulated  plus  one.  Each  row  of  ICL(I,J) 
represents  a discrete  point  in  time.  The  second  column  of 


ICL(I,J)  contains  a pointer  to  the  second  array,  the  event 
list  array  LIST(K).  The  first  column  of  ICL(I,J)  gives  the 
total  number  of  events  in  the  event  list  for  time  1.  The 
remaining  time  I events  follow  in  the  next  ICL(I,1)-1 
elements  in  LIST(K).  Time  is  incremented  as  discussed  above 
and  a modulo  function  is  used  to  determine  the  value  for  I 
from  TIME.  Thus,  the  simulation  time  at  which  an  event  is 
to  occur  is  described  implicitly  by  the  location  of  the 
pointer  in  the  ICL(I,J)  array.  The  use  of  these  two  arrays 
in  this  fashion  allows  a very  long  simulation  over  many 
increments  of  simulation  time  without  requiring  a large 
storage  array  to  maintain  timing  information. 

With  this  background,  a general  flow  through  the 
program,  can  be  outlined.  The  first  step  is  to  increment 
TIME.  The  clocking  array,  ICL(I,J),  is  then  accessed.  If 
there  are  events  scheduled  for  this  time,  they  are  executed 
as  listed  in  the  events  list  array.  Each  fanout  element  of 
each  element  in  the  event  list  is  in  turn  simulated  with  the 
appropriate  element  model.  An  element  status  and  descrip- 
tion array,  IDESOU(I,J),  is  used  to  maintain  a cell  of 
information  about  each  element  in  the  circuit.  The  informa- 
tion included  in  this  array  is  a list  of  fanout  elements  for 
each  element  of  the  circuit,  the  simulation  time  when  the 
last  output  change  occurred,  the  old  value  of  the  element 
output,  the  present  value  of  the  element  output,  the  simula- 
tion scheduled,  and  the  future  value  of  the  element  output. 
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With  this  information  the  element  model  determines  a new 
element  output  and  calculates  the  simulation  time  of  a 
future  element  output  change.  If  the  event  will  cause  a 
future  state  change  to  occur  in  the  fanout  element,  the 
future  event  is  properly  scheduled  into  the  event  list. 

When  all  events  have  been  processed,  time  is  again  incre- 
mented and  the  process  continues. 

It  is  interesting  to  note  that  if  the  cyclic  clock 
discussed  above  is  based  on  maximum  element  delay  time,  old 
information  in  the  event  list  will  automatically  be  replaced 
with  new  information.  This  reduces  the  required  length  of 
the  event  list.  For  example,  assxime  that  the  maximum  delay 
of  any  element  in  a circuit  is  10  units  of  time.  Define 
LIST (I)  to  have  I = 1000.  With  a clock  length  of  10  units, 
then  each  clock  unit  of  time  could  point  to,  on  an  average 
of  100  locations  in  LIST(IOOO)  or  100  events  could  be 
scheduled  per  increment  of  time.  Based  on  the  work  by 
Ulrich  (1969),  one  per  cent  of  the  circuit  elements  is  equal 
to  the  100  events  or  a circuit  comprising  of  10,000  elements 
could  be  simulated  by  this  example  without  increasing  the 
length  of  LIST (I) . 

It  has  already  been  pointed  out  that  an  element 
model  must  be  writtin  in  order  to  not  only  determine  the 
new  value  of  an  element,  but  also  to  determine  when  the 
change  is  to  occur.  The  information  required  to  make  these 


J 


calculations  must  include  the  simulation  time  when  the  last 
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output  change  occurred  and  the  old  value  of  the  output.  The 
model  must  also  know  if  the  element  has  already  been 
scheduled.  This  can  be  done  either  by  checking  the  event 
list  or  by  adding  this  information  to  the  above  and  storing 
it  with  the  element  description  cell  as  is  the  other  in- 
formation. Using  a cell  of  information,  of  course,  is  much 
more  desirable  than  searching  lists. 

Also  included  in  the  cell  of  information  is  the 
fanout  elements  of  the  element  concerned.  This  information 
is  used  to  determine  a new  set  of  active  elements.  In  order 
to  get  all  parameters  into  the  scheme,  internally  the 
program  generates  fictitious  elements  to  represent  each  of 
the  external  circuit  inputs.  The  program  also  expands 
multiple  output  elements  to  look  as  though  they  were  single 
output  elements.  Thus,  the  timing  scheme  does  not  need  a 
separate  mode  to  handle  multiple  output  elements.  All 
elements  and  the  external  circuit  inputs  are  all  handled  in 
the  same  fashion,  and  therefore  minimize  the  simulation 
control  effort. 


CHAPTER  4 


EVESIM  PROGRAM  DESCRIPTION 

EVESIM  has  been  written  in  FORTRAN  IV  for  use  on  the 
CDC  6400  computer  system.  EVESIM  is  best  described  as  being 
a time-driven  event  simulator  of  digital  logic  circuits. 

It  is  a fundamental  mode  simulator  in  that  elements  are 
modeled  in  terms  of  rising  and  falling  delay  times.  The 
program  has  been  divided  into  five  major  subroutines  as 
indicated  in  Fig.  1.  There  are  several  subroutines  that 
perform  functions  in  support  of  these  five  and  will  be  dis- 
cussed individually  as  their  time  of  use  is  identified.  As 
can  be  seen  in  Fig.  1,  the  first  three  subroutines  are  pre- 
processing routines.  Subroutine  FUNSIM  controls  the  actual 
simulation  of  the  logic  circuit  and  subroutine  OUTPUT 
controls,  as  the  name  implies,  data  output  for  the  program. 
The  program  has  debugging  aids  programmed  into  it  to  assist 
the  user  in  running  the  program  and  for  use  in  future  up- 
dating of  the  program. 

Preprocessing  Subroutines 

The  preprocessing  requirements  of  EVESIM  have  been 
broken  into  three  subroutines:  READIN,  INCOND,  and  INTIAL. 
Subroutine  READIN  reads  from  cards  the  circuit  description 
array  of  the  circuit.  Subroutine  INCOND  reads  from  cards 
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initial  circuit  values  and  initial  external  input  values, 
and  subroutine  INTIAL  initializes  the  circuit  in  preparation 
for  beginning  the  simulation. 

Subroutine  READIN 

Figure  2 describes  the  flow  of  functions  performed 
by  subroutine  READIN.  The  date  of  run,  the  circuit  number 
(LCKT) , the  total  number  of  simulation  times  to  be  executed 
(LAST) , and  the  mode  of  debug  (LBU)  are  read  from  the  first 
card.  Since  error  messages  can  be  generated  early  in  the 
program,  the  error  message  heading  is  printed  next.  The 
total  number  of  elements  in  the  circuit  is  read  in  as  NGATE, 
and  the  total  of  external  inputs  is  read  in  as  NEXIN. 

The  element  descriptions  are  read  in  using  the  same 
format  as  SCIRTSS.  The  circuit  elements  are  numbered  1 
through  NGATE  and  are  assigned  a 4-letter  type  designator 
such  as  JKMl  for  a J-K  master  slave  flip-flop.  The  input 
connections  to  the  element  are  then  listed  in  order.  The 
connection  is  indicated  by  the  number  of  the  element  whose 
output  is  connected  to  the  input  of  the  present  element. 

The  element  delay  information  follows  in  terms  of  simulation 
times  both  for  rising  and  falling  logic  values.  All  of  this 
information  is  stored  in  an  array  called  IDESCR(I,J).  Index 
I will  be  equal  to  NGATE,  the  number  of  elements,  and  Index 
J can  have  values  1 through  19  as  indicated  in  Table  1.  The 
element  type  designator  is  compared  against  the  available  ( 


Fig.  2.  Subroutine  READIN 
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Table  1.  Definition  of  J for  IDESCR(I,J). 


J Information  stored 


1 

2 

3 

4 

5 

6 

7 

8 

9-18 

19 


Element  type 

Element  type  number  used  for  internal  control 
Pointer  to  IDESOU (lOD, JO) 

First  value  of  delay  time  low  to  high 
First  value  of  delay  time  high  to  low 
Second  value  of  delay  time  low  to  high 
Second  value  of  delay  time  high  to  low 
Circuit  output  indicator 
Number  of  element  connected  to  input 
Status  of  element 
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element  models  listed  in  LTYAL(J)  and  the  appropriate  type 
number  assigned.  A fanout  description  of  the  circuit  is 
then  formed.  This  fanout  description,  based  on  the  fanout 
of  each  element,  is  in  effect  the  destination  set  discussed 
in  Chapter  2.  The  fanout  description  array  is  called 
IDESOU ( lOD, JO) . lOD  is  a function  of  the  number  of  elements 
in  the  circuit  and  the  number  of  outputs  of  each  element. 

An  element  that  has  more  than  one  output  is  represented  in 
IDESOU (lOD, JO)  once  for  each  output.  For  example,  each 
flip-flop  will  appear  twice,  once  for  the  Q output  and  once 
for  the  Q output.  The  JO  index  for  each  element  output  is 
defined  in  Table  2. 

Table  2.  Definition  of  JO  for  IDESOU (lOD, JO) . 

Information  stored 
Pointer  to  IDESCR(I,J) 

Element  type  number  used  for  internal  control 
Simulation  time  when  last  output  change  occurred 
Old  value  of  the  output 
Present  value  of  the  output 
Simulation  time  of  future  change 
Value  of  output  at  future  change 
Element  numbers  connected  to  the  output 
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Subroutine  INCOND 

Though  normally  all  required  data  for  a program 
would  be  read  into  the  computer  by  a single  subroutine,  this 
task  has  been  split  in  EVESIM.  This  allows  greater  flexi- 
bility in  input  data  formats.  Subroutine  READIN  is  oriented 
around  the  SCRITSS  input  format  and  may  not  be  suitable  to 
another  user's  needs;  therefore,  doing  preprocessing  in 
smaller  pieces  allows  easier  replacement  of  sections  to  meet 
a particular  user's  needs. 

Subroutine  INCOND  (Fig.  3)  reads  from  cards  the 
total  number  of  flip-flops  (LFF) , the  total  number  of  other 
memory  type  devices  (LOM) , two  minimum  values  of  clock 
pulse  widths  (NCLWID,  NCLW2) , and  two  minimum  values  of 
preset  or  clear  pulse  widths  (NPCWID,  NPCW2) . It  reads  the 
output  value  of  all  memory  devices  and  stores  that  informa- 
tion in  IDESOU (lOD, JOS) . The  final  data  read  by  INCOND  is 
the  initial  values  of  all  external  inputs. 

Subroutine  INTIAL 

The  most  important  part  of  the  preprocessing  is 
establishment  of  the  initial  values  for  all  elements  of  the 
circuit.  It  is  assumed  that  the  external  inputs  to  the 
circuit  have  remained  stable  for  a long  period  of  time,  in 
fact,  long  enough  for  any  previous  activity  to  have 
propagated  through  the  circuit.  This  situation  certainly 
cannot  be  achieved  if  an  unstable  condition  can  result  from 
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network  feedback  internal  to  the  circuit.  Such  an  unstable 
' condition  may  or  may  not  be  designed  into  the  circuit.  Any 

unstable  condition  is  recognized  by  the  program  and  the  user 
provided  with  sufficient  data  to  diagnose  the  situation. 

Subroutine  INTIAL  (Fig.  4)  begins  by  reading  the 
maximum  number  of  iterations  (LIMIT)  or  passes  through  the 
circuit  that  can  be  made  in  order  to  clear  all  unknowns  and 
establish  a stable  condition  in  the  circuit.  Failure  to 
establish  stability  will  result  in  error  messages  being 
printed  and  program  execution  stopped,  unless  of  course,  the 
initialization  override  parameter,  LOVER,  is  set  to  one. 

New  outputs  for  each  element  are  calculated.  The  unknown 
condition  X is  not  allowed  to  propagate  through  the  circuit, 
and  therefore  provides  a means  of  identifying  indeterminate 
situations  that  exist  in  the  circuit.  These  situations  may 
or  may  not  be  desired.  If  they  are  not  desired,  then  they 
have  been  identified  to  the  user.  Iterations  of  new  output 
calculations  are  continued  until  no  more  changes  occur  or 
the  limit  of  attempts  is  reached  as  outlined  in  the  flow 
chart  in  Fig.  4. 

Simulation  Control 

Efficiency  requirements  in  the  preprocessing  phase 
of  simulation  is  easily  put  aside  when  considering  that  the 
preprocessing  is  done  only  one  time.  Even  the  program 
efficiency  of  an  element  model  does  not  have  to  be  great  if 
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maximum  effort  is  exerted  to  call  the  model  only  when 
necessary.  Certainly  any  simulator  that  calculates  new 
element  outputs  for  each  element  during  each  increment  of 
time  requires  a very  efficient  element  model.  The  process 
of  searching  long  lists  as  proposed  by  pure  event  directed 
simulations  is  also  not  very  efficient.  EVESIM  does  not 
require  major  list  searching  routines  or  efficient  element 
modeling.  Great  flexibility  and  efficiency  have  been  gained 
in  EVESIM  by  not  involving  the  element  model  in  the  control 
and  driving  mechanism  of  the  simulator. 

Subroutine  FUNSIM  (Fig.  5)  is  the  control  and 
driving  mechanism  for  EVESIM.  The  event  list,  discussed  in 
Chapter  3,  is  defined  for  EVESIM  to  be  LIST(K).  K was 
chosen  to  be  equal  to  1024  for  purposes  of  easier  manipula- 
tion of  the  events  list;  however,  K is  by  no  means  limited 
to  this  number.  The  internal  clock  array  for  EVESIM  is 
defined  as  ICL(I,J)  and  simulation  time  is  defined  to  be 
equal  to  NOW. 

The  first  step  in  subroutine  FUNSIM  is  to  determine 
the  length  of  the  internal  clock,  LENCL,  which  is  equivalent 
to  the  longest  delay  time  that  has  been  assigned  to  any  one 
element  plus  one.  Then,  the  initial  value  of  the  elements' 
outputs  is  written  to  file  lOU,  and  the  first  event  card  is 
read.  An  event  card  contains  the  event  time  for  a change  in 
value  of  external  circuit  inputs,  the  external  input  number, 
and  the  new  value.  A card  can  have  only  one  time  value  but 
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may  contain  several  changes  for  that  time.  Additional  cards 

for  the  same  time  may  be  used  if  necessary.  The  present 

time  (LTIME)  of  the  clock  is  determined  by  a modulo 

function,  LTIME  = MODI (NOW, LENCL) , where  MODI  is  defined 

as  a function: 

FUNCTION  MODI  (M,N) 

LT  = M-(m:N)xN 
IF(LT.EQ.O)LT  = N 
MODI  = LT 

The  number  of  events  scheduled  for  a particular  time  is 
defined  as  NU.  NU  being  negative  in  value  indicates  that 
the  events  have  been  processed  for  that  time;  otherwise,  the 
scheduled  events  are  allowed  to  occur.  External  circuit 
input  changes  are  processed  and  then  internal  element 
changes  are  processed.  During  the  process  of  scheduling  the 
effects  of  element  changes,  each  fanout  element  of  the 
element  that  changed  is  processed  by  the  appropriate  element 
model  and  future  events  scheduled  as  needed.  Time  is  incre- 
mented and  the  process  continues  until  the  time  limit  is 
reached . 


Postprocessing  Subroutines 
Postprocessing  of  the  simulation  merely  provides  a 
usable  output.  The  process  has  been  divided  into  two  sub- 
routines, subroutine  OUTPUT  and  subroutine  PLOTl.  Sub- 
routine OUTPUT  is  used  for  control  purposes  while  subroutine 
PLOTl  does  the  actual  printing. 
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The  ext(irnal  circuit  inputs  and  the  element  outputs 
to  be  plotted  are  read  from  a card.  Up  to  12  outputs  can 
be  re<iuested  pc^r  card.  As  indicated  in  Fiq.  6,  the  output 
r(‘i|u<‘stf:  arf'  first  sorted  and  subroutine  IM.OTl  is  called  to 
do  the  I ) r i n t i luj . 

Subroutine  PLOTl  (Fig.  7)  reads  a record  from  the 
output  file  lOU.  It  then  fills  an  output  buffer  with  the 
appropriate  output  value  and  prints  the  contents  of  the 
buffer.  Another  record  is  read  from  lOU  and  the  process 
continues  until  the  file  is  empty.  It  should  be  obvious 
that  the  printing  routine  PLOTl  has  been  separated  from  the 
control  in  order  that  changes  to  the  actual  output  format 
can  be  easily  accomplished  with  minor  change  only  to  PLOTl 
and  not  the  whole  output  process. 

Supporting  Functions 

Discussion  to  this  point  has  only  referred  to  such 
things  as  modeling  routines,  scheduling  and  unscheduling 
events,  and  error  messages.  Though  subroutines  SCHED, 
UNSCHE,  ERROR,  and  lODEST  are  normally  referenced  only  by 
the  element  model  routines,  they  are  discussed  prior  to  the 
modeling  routines  (Chapter  4)  since  they  are  general  in 
nature  and  apply  to  all  the  models. 

Subroutine  SCHED 

The  purpose  of  subroutine  SCHED  (Fig.  8)  is  to 
schedule  an  event  to  occur  at  the  time  (NSCTI)  specified  by 


an  olomcnt  model.  It  is  a very  simple  process  and  involves 
dctermininq  whether  or  not  a block  has  been  reserved  in 
LIST(I)  for  time  NSCTI.  If  a block  exists,  it  is  opened  and 
the  event  inserted.  If  a block  does  not  exist,  one  is 
created. 

Subroutine  UNSCHE 

There  are  certain  instances  where  an  event  that  had 
been  previously  scheduled  needs  to  be  unscheduled.  To 
accomplish  this,  subroutine  UNSCHE  (Fig.  9)  is  used.  The 
event  to  be  removed  is  identified  by  the  element  model,  the 
pointer  (LPOIN)  to  the  beginning  of  the  block  in  LIST (I)  is 
then  established.  The  block  is  searched  for  the  event  and 
it  is  removed.  IDESOU (lOD, J06)  is  updated  to  reflect  that 
no  event  is  to  occur  and  control  is  returned  to  the  element 
model  subroutine. 

Subroutine  ERROR 

Subroutine  ERROR  (Fig.  10)  is  used  throughout  the 
program  to  process  error  messages  and  to  stop  program  execu- 
tion in  the  case  of  fatal  errors.  Errors  are  processed  by 
assignment  of  a number.  It  is  noted  by  Fig.  10  that  error 
messages  greater  than  199  are  fatal  and  cause  execution  to 
cease  by  using  the  library  function  CALL  EXIT.  After 
processing  of  non-fatal  errors,  control  is  returned  to  the 
calling  location.  A list  of  error  messages  that  can  occur 
is  contained  in  Chapter  7. 


luLroviCii 

UNSvm: 


Ifi»SOU(lODl,J06>  - 0 


9.  Subroutine  UNSCHE 
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Subroutine  lODEST 

Subroutine  lODEST  establishes  the  index  lOD  of 
IDESOU (ion, JO)  for  the  current  element  being  processed 
during  the  simulation.  It  establishes  lODl  as  the  element 
being  proeessed,  I0D2  as  the  first  input  to  the  element, 

I0D3  as  the  second  input,  and  on  to  as  many  indices  as  there 
are  inputs  to  the  element.  It  is  merely  a "look-up"  routine 
intended  to  minimize  cross-referencing  efforts  between 


IDESOU  and  IDESCR  arrays. 


CHAPTER  5 


EVESIM  ELEMENT  MODELS 

Many  approaches  exist  to  model  logic  elements. 
Chapter  2 established  that  an  element  level  approach  was  the 
most  appropriate  since  modeling  elements  with  gates  may  not 
present  the  same  operating  characteristic  as  the  element 
itself.  It  was  also  established  that  an  assignable  or 
precise  timing  capability  was  most  desirable  in  order  to 
achieve  great  accuracy.  In  fact,  being  able  to  assign 
delays  in  terms  of  both  rising  signal  levels  and  falling 
levels  provides  the  most  accurate  model.  EVESIM  has  been 
designed  such  that  delays  can  be  assigned  in  terms  of  rising 
and  falling  signal  levels.  The  problem  of  modeling  for 
EVESIM  is  a matter  of  assigning  the  correct  delay  times  to 
the  element  and  being  able  to  calculate  from  these  delays 
when  an  element's  outputs  are  to  change.  The  models  have 
been  split  into  two  distinct  types.  The  first  type  is  for 
single  output  non-memory  type  devices  normally  called 
combinational  logic  gates.  The  second  type  is  for  all  other 
devices  which  were  defined  earlier  as  elements. 

Element  Delay  Time 

EVESIM  uses  precise  delay  time  calculations. 
Obviously,  either  efficiency  or  accuracy  would  quickly  be 
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lost  if  element  delay  times  were  randomly  picked  or  an 
attempt  was  made  to  run  the  simulation  in  real  time.  It  is 
therefore  necessary  to  establish  certain  ground  rules  for 
determining  the  element  delay  times. 

It  is  important  that  the  delay  times  for  the 
elements  have  some  common  parameter;  otherwise,  the  accuracy 
of  the  simulation  is  lost.  If  all  element  delays  were  made 
to  be  of  the  form:  D = nAt,  where  D is  delay  time  in  real 

time,  n is  a positive  integer,  and  At  is  a positive  constant, 

then  all  delays  can  be  represented  by  a positive  integer  n. 

This  has  forced  all  element  delays  in  a given  circuit  to  be 

relative  to  each  other;  and  therefore,  the  accuracy  of  the 
simulation  has  been  preserved  while  not  losing  efficiency. 

An  EVESIM  user  would  refer  to  a manufacturer's  logic  data 
book  and  determine  approximate  values  for  t „ and  t for 

IT  Ltii  r^HXj 

all  the  elements  in  the  circuit  to  be  simulated.  It  should 
be  noted  here  that  the  choice  of  delay  times  is  not  totally 
random  and  requires  the  user  to  have  some  intuitive  knowl- 
edge of  the  expected  operation  of  the  circuit  being  simu- 
lated. Then,  from  this  information,  the  user  should  deter- 
mine a value  for  At  which  by  definition  is  nothing  more  than 
the  greatest  common  divisor  of  the  element  delays.  All 
element  delay  times  (n)  can  then  be  determined  by:  “ 


D : At . 
e 


As  an  example,  consider  a small  circuit  containing 
two  types  of  elements  with  delays  of  = 10  nsec. 


45 


^PHLl  ~ ^ nsec,  tpj^p^2  “ nsec,  and  tp^^2  “ ^ nsec.  The 
greatest  common  divisor  of  these  times  is  2 nsecs  which  is 
equal  to  At.  The  values  for  n are  then  computed: 


ILH  = ^ 

2=5 

"^IHL  ^ ■ 

2 

2HL  " ' 

: 2 = 6 

'^2HL  ""  ® ^ 

2 

One  simulation  time  is  now  equivalent  to  2 nsec.  The 

maximum  length  of  the  internal  control  clock  has  been 

determined  and  is  equal  to  the  maximum  n which  in  this  case 

e 

is  6.  It  should  be  noted  that  another  set  of  delays  will 
determine  a different  At.  Thus,  it  is  desirable  to  maximize 
At  such  that  the  maximum  element  delay  time  n^  is  kept  as 
small  as  possible  in  order  to  economize  the  execution  of  a 
simulation. 


Subroutine  COMLOG 

When  an  event  has  occurred,  an  element  output  has 
changed.  The  elements  connected  to  that  changed  output  are 
each,  in  turn,  processed  through  an  element  model.  Sub- 
routine COMLOG  (Fig.  11)  serves  as  the  model  for  all  com- 
bintational  logic  gates.  The  first  action  taken  is  a calcu- 
lation of  the  output.  Then,  a comparison  is  made  to  the  old 
output.  If  a change  has  occurred,  the  new  output  is 
appropriately  scheduled  as  an  event  to  occur  in  the  future. 

There  are  many  techniques  available  that  can  be  used 
to  determine  the  logic  values.  Boolean  algebra  can  be  used. 
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The  available  logic  functions  of  an  operating  system  provide 
a means  of  calculation.  Table  look-up  techniques  can  be 
used  or  arithmetic  expressions  can  be  used  for  the  calcula- 
tions. The  technique  used  for  the  calculation  is  not  really 
relative  to  the  modeling  of  the  element  as  far  as  modeling 
the  delay  is  concerned.  Therefore,  whatever  is  easiest  for 
the  programmer's  system  is  appropriate.  EVESIM  uses  both 
table  look-up  techniques  and  the  available  logic  functions. 

The  important  part  of  the  modeling  scheme  is  that  an 

event  that  is  to  occur  is  properly  scheduled  in  the  future. 

This  scheduling  is  based  on  the  delay  times  as  calculated 

above.  It  is  important  that  noise  spikes  are  passed  or 

rejected  as  required.  Whether  a gate  will  pass  or  reject  a 

spike  is  a function  of  the  gate  propagation  delay  times, 

tpLH  srid  The  modeling  of  this  is  easily  accomplished 

with  an  event  type  of  simulation  program.  For  example,  let 

a NAND  gate  have  values  t„^„  = 11  and  t„„^  = 7.  Obviously, 

immunity  to  1 - to  - 0 transitions  is  greater  than  0 - to  - 

1 transitions  since  is  greater  than  Let  the 

output  of  the  NAND  gate  at  time  0 be  equal  to  0.  Assume 

that  at  time  0,  an  input  to  the  gate  has  changed  such  that 

the  output  should  change  to  1.  An  event  should  be  scheduled 

for  time  0 + tn^„  which  is  equal  to  time  11.  Following  the 

scheme  of  the  flow  chart  in  Fig.  11,  any  attempt  to  schedule 

a t^uT  transition  prior  to  time  11  would  negate  the  event 
r HL 


scheduled  at  time  11.  The  earliest  that  a valid  t 
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transition  would  reqister  would  be  at  time  = 11  + t^„-  = 18. 
Therefore,  the  narrowest  pulse  that  could  be  passed  through 
the  NAND  gate  would  be  18  - 11  = 7,  which  is  what  is 
expected.  A similar  analysis  can  be  done  for  the  1 - to  - 0 
noise  immunity.  The  result  is  a minimum  pulse  of  11,  again 
which  is  what  is  expected.  Making  use  of  the  rising  and 
falling  delay  times  in  conjunction  with  proper  scheduling 
and  unscheduling  of  events  has  effectively  modeled  the 
combinational  logic.  The  modeling  of  multiple  output 
devices  is  accomplished  in  a similar  fashion.  The  differ- 
ence, as  seen  in  the  next  section,  is  in  the  complexity  of 
determining  when  an  event  should  be  scheduled. 

Subroutine  JKMl 

The  complexity  of  determining  when  an  event  should 
occur  for  elements  such  as  a JK  master-slave  flip-flop 
results  from  the  fact  that  the  flip-flop  has  several  inputs, 
it  has  memory  capability,  and  has  more  than  one  output. 

There  are  also  other  parameters  for  consideration  such  as 
minimum  pulse  width  requirements  to  the  clock  input  and  the 
preset  and  clear  inputs,  and  the  J and  K inputs  remaining 
stable  while  the  clock  is  high.  Since  the  data  depicting 
these  parameters  is  available  in  the  data  book,  they  are 
also  applied  against  the  greatest  common  divisor  and  an 
integer  value  assigned  to  represent  the  parameters  for 


modeling  purposes. 
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The  JK  flip-flop  model,  subroutine  KJMl  (Fig.  12), 
is  broken  into  four  paths  dependent  on  the  three  inputs  that 
can  possibly  change  the  outputs.  Since  preset  and  clear 
override  the  clock,  they  first  are  tested  for  an  unknown 
condition.  An  unknown  on  preset  or  clear  will  automatically 
force  an  unknown  condition  on  the  output.  If  neither  are 
unknown,  further  checks  are  made  and  the  appropriate 
branching  takes  place  as  indicated  in  Fig.  12.  If  the  J or 
K input  to  the  flip-flop  was  the  input  change  that  caused 
the  model  to  be  called,  control  is  returned  to  the  calling 
subroutine  without  any  further  action.  It  should  also  be 
noted  that  a leading  edge  clock  change  requires  no  action 
for  a master-slave  flip-flop  and  control  is  returned  to  the 
calling  routine. 

Assuming  that  a trailing  edge  clock  caused  the  call 
of  subroutine  JKMl,  the  clock  pulse  is  first  checked  to 
insure  that  it  meets  minimum  width  requirements.  If  not, 
an  error  message  (error  1)  is  processed  and  the  outputs 
scheduled  as  unknowns.  If  the  clock  is  valid  and  preset  and 
clear  are  both  high,  new  outputs  are  calculated  and 
scheduled  to  occur  as  events  in  the  future. 

A single  preset  change  or  a clear  change  are  similar 
in  nature  and  relatively  simple  to  program.  The  processing 
of  these  changes  are  self  explanatory  and  are  indicated  on 
the  flow  chart  on  pages  52  and  54 . The  complexity  is 
greatly  increased  if  both  preset  and  clear  change  at  the 
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same  time.  If  they  were  both  low,  the  outputs  are  both 
high.  If  both  preset  and  clear  go  high  at  the  same  time, 
it  can  not  be  determined  what  the  outputs  will  be;  thus,  the 
outputs  are  changed  to  the  unknown  condition.  If  preset  and 
clear  are  opposite  in  value  and  they  both  change,  a spike 
can  be  generated  on  the  output  of  the  flip-flop.  This  is 
designated  by  error  message  6 and  the  outputs  are  tempo- 
rarily changed  to  an  unknown  condition  until  the  new  values 
of  preset  and  clear  cause  the  outputs  to  stabilize. 

On  the  surface,  modeling  of  the  JK  master-slave 
flip-flop  appears  to  be  a lengthy  process.  However,  by 
having  full  understanding  of  how  the  element  functions,  the 
model  can  be  broken  into  small  processing  sections  each 
representing  an  element  reaction  to  a particular  stimulus, 
thus  providing  an  easily  understood  and  logical  approach  to 
solving  the  modeling  problem  for  a JK  flip-flop  or  any 
multiple  input-multiple  output  element. 
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CHAPTER  6 

DEMONSTRATION  OF  THE  USE  OF  EVESIM 

A description  of  EVESIM  would  not  be  complete  with- 
out a demonstration  of  some  circuits  being  simulated.  The 
first  example  is  strictly  a JK  master-slave  flip-flop  run 
against  the  simulator  all  by  itself.  This  was  chosen  in 
order  to  demonstrate  the  functioning  of  the  flip-flop  model 
discussed  in  Chapter  4.  The  remaining  two  circuits  were 
chosen  to  demonstrate  the  interaction  of  the  models.  It 
should  be  noted  that  both  have  critical  feedback  paths  that 
restrict  the  frequency  of  the  clock,  thus  surely  demon- 
strating the  need  for  a level  mode  simulation.  Data  used 
for  the  simulations  were  taken  from  the  Texas  Instruments 
TTL  Data  Book  (Components  Group,  Engineering  Staff,  1973). 

Demonstration  of  the  JK  Master-Slave 
Flip-Flop 

The  flip-flop  operation  has  been  demonstrated  by  six 
different  runs  of  the  simulator.  The  six  runs  were  designed 
in  order  that  all  possible  paths  through  the  flip-flop  model 
would  be  exercised.  Figure  13  identifies  the  input  names  to 
the  flip-flop  and  Table  3 contains  the  data  used  for  the  six 
runs.  It  should  be  noted  that  the  output  of  the  flip-flop 
is  given  by  the  Q output  being  equal  to  element  number  1 and 
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Table  3.  JK  master-slave  flip-flop  data. 

1.  Greatest  common  divisor  = At  = 5 nanoseconds 

2.  Clocked  Input  Data 

t„^ „ = 15  nanoseconds  n = 15  t 5 = 3 Simulation  times 
PLH 

t = 25  nanoseconds  n = 25  v 5 = 5 Simulation  times 

IT  11  Jj 

Minimum  allowed  clock  pulse  width  = 20  nanoseconds 
NPCWID  = 20  T 5 = 4 Simulation  Times 


Preset 

and 

Clear  Input 

Data 

'^PLH 

20 

nanoseconds 

n = 

20  r 

5=4 

Simulation 

times 

^PHL 

30 

nanoseconds 

n = 

30  T 

5 = 6 

Simulation 

times 

Minimum  allowed  preset  or  clear  input  pulse  width  = 
25  nanoseconds 

NPPWID  = 25  5 = 5 Simulation  Times 
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the  Q output  by  element  number  2001.  The  2000  indicates  the 
second  output  of  element  number  1 or  2001. 

Figure  14  contains  the  first  run.  The  J input  high 
and  the  K input  low  at  time  2 causes  no  change  at  the 
outputs  as  is  expected.  Similarly,  with  both  J and  K low, 
the  clock  pulse  at  time  14  causes  no  change  at  the  output. 
The  clock  pulse  at  time  26  causes  the  output  to  change  since 
the  K input  is  high.  It  should  also  be  pointed  out  that  the 
model  at  this  time  is  demonstrating  the  fact  that  a JK 
master-slave  flip-flop  requires  a set  up  time  equal  to  zero 
and  a hold  time  equal  to  zero.  In  other  words,  the  J and  K 


inputs  must  be  stable  while  the  clock  is  high  in  order  to 
guarantee  proper  operation  of  the  master-slave  flip-flop. 
With  both  J and  K inputs  high,  the  clock  pulse  at  time  38 
causes  the  output  to  toggle. 

Figure  15  is  the  output  of  the  flip-flop  for  run 
number  2.  Run  2 is  intended  to  demonstrate  the  operation  of 
the  clear  and  preset  inputs.  Notice  that  at  time  1 the 
clear  input  goes  low,  overriding  the  clock  pulse  and  causing 
an  output  change.  The  preset  is  similarly  demonstrated  at 
time  11.  Notice  that  at  time  18  both  preset  and  clear  are 
low  causing  both  outputs  to  be  high.  At  time  23  the  clear 
input  returns  to  the  one  level  and  the  outputs  align  them- 
selves according  to  clear  remaining  low. 

The  remaining  four  runs  made  of  the  JK  flip-flop 


demonstrate  the  possible  error  conditions  that  can  cause  an 
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TIME 


CIRCUIT  NUMBER  » 3 

ELEMENT  NUMBER 


TIME 


2001 

1 

1001 

1002 

1003 

1008 

1005 

0 

0 

1 

0 

1 

1 

1 

1 

0 

1 

0 

1 

0 

1 

1 

0 

1 

1 

2 

0 

1 

0 

1 

1 

0 

1 

2 

3 

0 

1 

1 

1 

1 

0 

1 

3 

0 

1 

1 

1 

1 

0 

1 

8 

1 

1 

1 

1 

1 

0 

1 

5 

6 

1 

1 

1 

1 

1 

1 

1 

6 

7 

1 

0 

0 

1 

1 

1 

1 

7 

8 

1 

0 

0 

1 

1 

1 

1 

8 

1 

0 

0 

1 

1 

1 

1 

9 

10 

1 

0 

0 

1 

1 

1 

1 

10 

1 1 

1 

0 

0 

1 

1 

1 

0 

11 

12 

1 

0 

0 

1 

1 

1 

0 

12 

U 

1 

0 

0 

1 

1 

1 

0 

13 

18 

1 

0 

0 

1 

1 

1 

0 

18 

15 

1 

1 

1 

1 

1 

1 

0 

15 

16 

1 

1 

1 

1 

1 

1 

1 

16 

1 7 

0 

1 

1 

1 

1 

1 

1 

17 

18 

0 

1 

1 

1 

1 

0 

0 

18 

19 

0 

1 

0 

1 

1 

0 

0 

19 

20 

0 

1 

0 

1 

1 

0 

c 

20 

21 

0 

1 

0 

1 

1 

0 

0 

21 

22 

1 

1 

0 

1 

1 

0 

0 

22 

23 

1 

1 

0 

1 

1 

0 

1 

23 

28 

1 

1 

0 

1 

1 

0 

1 

28 

2 5 

1 

1 

0 

1 

1 

0 

1 

25 

26 

1 

1 

0 

1 

1 

0 

1 

26 

27 

1 

1 
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1 

0 

1 

27 
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0 

1 

28 

29 
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1 

1 

29 
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1 

1 
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1 
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1 

1 

1 

31 

32 

1 

0 
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1 

1 
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1 

0 

0 

1 

1 

1 

1 
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38 

1 

0 

0 

1 

1 
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unknown  condition  to  appear  on  the  output.  The  error 
messa<jes  are  clearly  printed  with  the  error  number,  time  of 
occurrence,  and  the  element  involved.  Figure  16  demonstrates 
two  detected  errors.  The  first  is  error  1,  which  is  a clock 
pulse  being  too  narrow.  This  error  occurred  at  time  2. 

Error  7,  though  not  necessarily  an  error  in  the  true  sense 
of  the  word  error,  points  out  that  the  outputs  are  unknown, 
and  since  the  J and  K inputs  are  of  the  same  value,  a new 
set  of  outputs  cannot  be  determined.  Notice  that  at  time 
23,  the  J and  K inputs  have  different  values,  and  therefore 
a new  output  set  can  be  determined.  The  worst  case  is 
assumed;  that  is,  that  for  the  Q output  the  unknown  was 
equal  to  1 and  had  to  change  to  0,  and  that  the  opposite 
case  exists  for  Q.  The  remaining  point  of  interest  on  Fig. 

16  is  that  at  time  35,  a clock  pulse  occurs  with  the  Q out- 
put low  and  the  K input  low,  and  because  of  this  time 
situation,  the  K input  changing  midway  through  the  clock 
pulse  has  no  effect  on  the  output  as  should  be  expected. 

At  time  5 on  Figure  17,  error  message  is  generated. 
This  error  results  from  ^he  fact  that  the  K input  to  the 
flip-flop  does  not  meet  the  hold  time  required,  or  the  K 
input  is  unstable  while  the  clock  is  high.  This  results  in 
the  unknown  condition  being  generated  on  the  outputs.  The 
next  interesting  point  is  the  spike  generated  on  the  Q out- 
put at  time  20.  This  is  a result  of  the  clear  input  having 
a low  value  and  the  preset  input  going  low  temporarily  and 
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then  returning  to  the  high  value.  Certainly,  this  operation 
can  be  expected.  Error  message  4,  clear  input  pulse  is  too 
narrow,  occurs  at  time  32  because  the  clear  input  was  not 
low  a sufficient  amount  of  time  after  the  preset  input  re- 
turned to  one.  Error  4 occurs  again  at  time  40. 

Figure  18  is  the  output  for  run  5.  There  are  3 
errors  demonstrated  by  this  run.  Error  3 occurs  when  both 
preset  and  clear  inputs  change  to  a high  value  at  the  same 
time.  Error  5 appears  when  a spike  occurs  on  the  preset 
input,  in  other  words,  preset  was  low  an  insufficient  amount 
of  time  to  cause  a predicted  change  on  the  output.  Error  7 
occurs  again  for  the  scime  reason  as  stated  earlier. 

Error  number  6 is  demonstrated  in  Fig.  19,  which  is 
run  number  6 of  the  flip-flop.  Error  6 is  generated  because 
preset  and  clear  are  changing  at  the  same  time  before  either 
preset  or  clear  was  able  to  cause  a change  to  the  output  of 
the  flip-flop.  Notice  that  this  situation  is  self- 
correcting  if  either  clear  or  preset  remains  low  for  a 
sufficient  amount  of  time. 

In  the  following  sections,  the  JK  master-slave  flip- 
flop  model  is  used  in  conjunction  with  the  combinational 
logic  model  to  simulate  realizable  circuits.  Having  fully 
established  the  operation  of  a single  element,  it  is  now 
time  to  demonstrate  the  simulation  of  a complete  circuit. 
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U ‘inonstration  of  an  Odd-Even 
Detector  Circuit 

Fiqure  20  is  a diaqram  of  an  odd-even  detector.  It 
is  a sequential  circuit  with  two  pulse  inputs,  A and  B,  and 
a sinqle  level  output  z.  Following  a pulse  on  line  B,  the 
output  z is  to  be  1 provided  there  has  been  an  even  number 
of  pulses  on  line  A since  the  previous  pulse  on  line  B. 
Otherwise,  a B pulse  will  reset  the  output  to  0.  The  output 
will  not  change  except  on  the  arrival  of  a B pulse.  The 
intent  of  the  example  is  to  determine  whether  the  circuit 
presented  is  or  is  not  the  best  design,  and  in  turn  to 
demonstrate  the  usefulness  of  a level  mode  simulator  like 
EVESIM. 


The  element  delay  parameters  used  for  the  first  run 
are  in  Table  4.  The  circuit  diagram  is  contained  in  Fig. 

20.  Under  the  conditions  of  Table  4,  the  circuit  functions 
as  it  should  and  the  output  of  the  simulator.  Fig.  21, 
verifies  that  it  does  function  correctly.  However,  the 
circuit  does  not  function  correctly  if  the  pulse  on  line  B 
is  made  wider  and  gate  number  4 is  made  faster.  This  situa- 
tion causes  flip-flop  2 to  change  value  very  fast;  thus, 
the  inputs  to  flip-flop  1 will  be  unstable  while  the  clock 
pulse  to  flip-flop  1 is  still  high.  The  parameters  for 
gate  4 are  changed  to  t„  „ = 10  nanoseconds  and  t = 5 

r Ijii  Hrl  Jj 

nanoseconds  and  pulse  B is  widened  to  45  nanoseconds.  With 


Table  4.  Odd-even  detector  circuit  parameters. 


Greatest  Common  Divisor  = At  = 5 nanoseconds 

Flip-Flop  Data 

a.  Clocked  input  data 

t = 25  nanoseconds  n = 25  v 5 = 5 Simulation 

Times 

t = 40  nanoseconds  n = 40  t 5 = 8 Simulation 

Times 

Minimum  allowed  clock  pulse  width  = 20  nanoseconds 
NPCWID  = 20  T 5 = 4 Simulation  Times 

b.  Reset  Input  Data 

t = 15  nanoseconds  n = 15  t 5 = 3 Simulation 

Times 

t = 25  nanoseconds  n = 25  ^ 5 = 5 Simulation 
P HJ-I  m • 

Times 

Minimum  allowed  reset  pulse  width  = 25  nanoseconds 
NPPWID  = 25  T 5 = 5 Simulation  Times 

Inverter  Data 

t„^ „ = 15  nanoseconds  n = 15  t 5 = 3 Simulation 

Times 


10  nanoseconds  n 


10  5 


2 Simulation 
Times 
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these  changes,  which  are  not  unreasonable,  the  circuit  no 
longer  functions  correctly  as  indicated  in  Fig.  22. 

Demonstration  of  an  Up-Down 
Gray  Code  Counter 

The  final  circuit  to  be  demonstrated  is  an  up-down 
self-correcting  gray  code  counter.  The  circuit  has  three 
input  signals,  a control  signal  X,  a control  signal  K,  and 
a clock  input.  When  X = 1,  the  counter  counts  up,  and  when 
X = 0,  the  counter  counts  down.  The  K input  controls  when 
the  counter  shall  be  allowed  to  count.  K = 0 does  not  allow 
counting,  and  K = 1 allows  counting  to  occur.  A circuit 
diagram  of  the  counter  is  presented  in  Fig.  23.  Typical 
values  for  delays  were  assigned  as  indicated  in  Table  5. 

Two  runs  were  made.  The  first  indicates  proper  functioning 
and  the  second  demonstrates  a fatal  user  error. 

Figure  24  contains  the  results  of  the  simulation  of 
the  counter  as  it  performs  all  its  required  functions.  The 
circuit  was  driven  with  a clock  pulse  equal  to  approximately 
8 mHz  with  a pulse  duration  of  36  nsec.  Attempts  to  drive 
the  circuit  with  a faster  clock  failed,  thus  indicating  that 
this  circuit's  maximum  courting  frequency  would  be  about 
8.0  mHz.  Assuming  worst  case  delays  and  adding  along  the 
longest  propagation  paths.  Gate  10  (22  nsecs)  + Gate  11 
(40  nsecs)  + Gate  7 (22  nsecs)  + Gate  8 (22  nsecs)  + Gate  12 
(25  nsecs)=  131  nsecs.  This  would  result  in  a predicted 
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Fig.  22.  Circuit  5 run  2 
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Table  5. 

Up-down  gray  code  counter 

circuit  parameters. 

Gate 

JK  Flip-flop 
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NSEC) 
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Fig,  24 


Circuit  4 run  1 


78 


79 


1!>2 

113 

15'. 

L55 

156 

15/ 

159 
tbC 
lol 
162 
163 
169 
lo5 
1 o6 
167 
166 
169 
1 7C 
1 71 
I 7? 
173 
I 79 

175 

176 

177 
L 76 

179 

180 
181 
162 
183 
[o'. 

195 
186 
18  7 
186 
169 

190 

191 

192 

193 
19'. 
1 <5 

196 

197 
196 
199 
2ClC 
201 
202 


0 

0 

0 

0 

0 

0 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

l_ 

I 

1 

1 

1 

1 

1 

X 

1 

I 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1 0 

1 0 

1 0 

I 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

I 0 

1 0 

I 0 

1 0 

1 0 

1 0 

. 0 

I 0 

I 0 

1 -J 

1 0 

1 0 

I 0 

I 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

1 0 

I 0 

1 0 

i 0 

1 u 

1 0 

1 0 

1 0 


1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

c 

c 

0 

u 

0 

c 

0 

0 

0 

0 

0 

c 

0 

c 

c 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

c 

0 

0 


1 

1 

c 

c 

0 

c 

0 

c 

0 

0 

0 

0 

0 

0 

0 

c 

0 

1 

1 

1 

1 

1 


c 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

c 

I 

1 

1 

1 

1 

1 

0 

0 

0 

c 

0 

0 

0 


152 

153 

159 

155 

156 

157 
1 56 
:59 

160 
161 
162 
163 
169 

165 

166 
lo7 
1 66 

169 

170 

171 

172 

173 
179 
175 
ITe 
177 
1 78 
] 79 
160 
Ibl 
182 
163 
169 
135 
166 
187 
166 
169 

190 

191 

192 

193 
199 

195 

196 

197 

198 

199 
2U0 
201 
202 


Fig.  24 . -^Continued  Circuit  4 run  1. 
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clock  frequency  of  7.6  mHz.  Certainly,  the  simulation  is 
valid. 
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Figure  25  is  provided  in  order  to  point  out  what 
happens  when  the  user  makes  an  input  control  error.  Notice 
that  the  user  has  requested  output  for  element  number  35. 
Referring  back  to  Fig.  23,  it  can  be  seen  that  no  such 
element  existed. 


CIRCUIT  NUMBER  * ^ 

MESSAGES************** 

ERROR  = 2C5  THIS  IS  A FATAL  ERROR 

LIST  OF  REQUESTED  ELEMENT  OUTPUTS 

1 2 3 A 5 35  1001  -0  -0  -0  -0  -0 


Fig.  25.  Circuit  4 run  2. 


CHAPTER  7 


USER'S  GUIDE  TO  EVESIM 

This  chapter  describes  the  use  of  EVESIM  as  it 
exists  at  The  University  of  Arizona.  It  outlines  the 
control  cards  required  and  the  input  data  format.  Also 
included  is  a detailed  list  of  the  existing  error  messages 
and  their  interpretation.  The  last  two  sections  of  the 
chapter  are  dedicated  to  explaining  the  debug  capabilities 
of  the  program  and  explaining  how  a user  would  write  his 
two  element  models. 


General  Input  Data 

There  are  two  phases  to  the  process  of  using 
EVESIM.  The  circuit  to  be  simulated  must  be  prepared  for 
the  simulation  and  the  required  data  must  be  punched  on 
cards.  Though  the  circuit  preparation  requirements  have 
been  discussed  throughout  the  previous  chapters,  it  is 
repeated  for  user  clarity. 

Circuit  Preparation 

EVESIM  must  be  provided  with  a description  of  the 
actual  circuit  to  be  simulated.  The  information  is  taken 
from  the  circuit  diagram.  This  requires  that  each  element 
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have  a unique  number  for  identification.  The  numbering  is 
begun  at  .1  and  continues  sequentially  through  n elements. 

An  EVESIM  user  must  then  refer  to  a manufacturer's 
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logic  data  book  and  determine  the  values  for  t__„  and  t„„^ 
for  all  the  elements  in  the  ::ircuit.  From  this  information 
the  user  should  determine  a value  for  At  which  by  definition 
is  the  greatest  common  divisor  of  the  element  delays.  All 
element  delay  times  are  then  determined  by 


where 

D is 
e 

the 

element 

delays  tp^^^  ^ 

or  tp^j^  as  appropriate. 

As  an 

example , 

consider 

two  types  of 

elements 

with  delays 

^PLHl 

= 10 

nsec 

' ^PHLl 

■ 6 nsec,  “ 12  nsec,  and 

^PHL2 

= 8 nsec. 

The  greatest  common 

divisor 

(At)  is  2 nsec 

The  values 

for 

n are  then  computed: 

n 

ILH 

10 

2=5 

’^IHL 

6 V 2 = 

3 

n 

2LH 

12  T 

2 = 6 

'^2HL 

II 

CN 

00 

4 

One  simulation  time  is  now  equivalent  to  2 nsec  and  all 
delay  parameters  have  been  defined  in  terms  of  simulation 
times.  It  should  be  noted  that  all  timing  information  now 
required  as  input  data  shall  be  in  this  same  format, 
adjusted  by  At. 

Card  Input  Format 

The  data  cards  are  presented  and  discussed  in  the 
sequence  that  they  should  appear  for  a proper  program 


execution.  The  normal  system  control  cards  required  are: 
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JOB  CARD,  CM70000 
ATTACH  (SIM,  ID  = MOEN) 

SIM 

7/8/9 

The  circuit  data  cards  follow  the  system  control 
cards.  Unless  otherwise  noted,  all  numbers  are  right 
justified  integers  in  5 column  fields. 

Card  1:  Title  Card.  The  first  10  columns  are  used  for 
the  month,  day,  and  year;  for  example,  Marl0,1975. 
The  circuit  number  is  right  justified  to  column  20. 
The  total  number  of  simulation  times  to  be  executed 
is  right  justified  to  column  25  and  the  debug  mode 
indicator  is  right  justified  to  column  30.  The  use 
of  debug  is  discussed  later. 

Card  2:  Circuit  Parameters.  NGATE  is  the  total  number 
of  elements  in  the  circuit  (col.  1-5).  NEXIN  is  the 
total  number  of  external  inputs  to  the  circuit  (col. 
6-10)  . 

Card  3:  Element  Description.  Each  circuit  element 

description  is  contained  on  a single  card.  That  is, 
each  element  in  the  circuit  will  have  a card  with 
the  following  format: 

Element  circuit  number  (col.  1-5)  . 

Element  type  (col.  7-10) . The  element  type  must  be 
one  of  the  following:  JKMl , NAN2,  NAN3,  NAN4, 
NAN8,  NOR2,  INVE,  XOR2 , OR2 , or  AND.  All 


r 
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flip-flops  are  to  be  listed  first.  Then  all 
other  memory  devices  followed  by  the  combina- 
tional logic  elements. 

Element  input  connections  to  this  element.  They  are 
listed  in  5-column  fields  beginning  with  columns 
16-20  up  to  column  45  and  then  beginning  with 
columns  56-60  up  to  column  75.  The  maximum 
number  of  inputs  is  10.  External  input  X is 
designated  by  1000  + X.  A constant  logical  1 is 
indicated  by  1000.  A constant  logical  0 is 
indicated  by  -1000.  The  second  output  of  an 
element  Y,  for  example  Q of  a flip-flop,  is 
indicated  by  2000  + Y.  The  third  output  would 
be  indicated  by  3000  + Y,  etc. 

Element  delay  information  appears  as  indicated  below 
(right  justified) : 

cloc)ced  or  normal  t„Tu  (col.  47-48) 
cloclced  or  normal  t^u^  (col.  49-50) 
preset,  clear  or  second  (col.  51-52) 

preset,  clear  or  second  t (col.  53-54) 

PHIj 

Card  4:  Memory  Element  Information. 

LFF  is  the  total  number  of  flip-flops  (col.  1-5) 

LOM  is  the  total  number  of  other  memory  type 
elements  (col.  6-10) 

NCLWID  is  the  first  "'inimvr'  c’.oc)c  pulse  width 


(col.  11-15) 
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NPCWID  is  the  first  minimum  preset  or  clear  pulse 
width  (col.  16-20) 

NCLW2  is  the  second  clock  pulse  width  (col.  21-25) 
NPCW2  is  the  second  minimum  preset  or  clear  pulse 
width  (col.  25-30) 

Card  5:  Initial  value  of  the  Q output  of  all  flip-flops, 
listed  in  order  as  they  appear  above,  beginning  with 
column  one.  Every  column  is  used.  A 0 represents 
logical  0.  A 1 represents  logical  1.  A 2 repre- 
sents an  unknown  condition. 

Card  6:  Same  as  card  5 except  for  all  other  memory 
devices,  if  any,  otherwise  skip  this  card. 

Card  7:  Initial  value  of  external  inputs  listed  in 
order,  beginning  with  the  clock  (1001).  Again, 
each  column  should  have  a 0,  1,  or  2. 

Card  8:  Initialization  parameters  LIMIT  is  the  limit  of 
the  number  of  attempts  that  will  be  made  to 
initialize.  LOVER  is  set  to  1 if  the  user  desires 
to  override  LIMIT,  otherwise  LOVER  is  left  blank. 

Card  9:  External  Input  Changes. 

lETIME  is  the  simulation  time  that  the  indicated 
external  inputs  are  to  change  value  (col.  1-5) 
INUV(15,2)  is  the  number  (1000  + K)  of  the  external 
input  that  is  to  change  and  the  value  (0,  1,  or 
2)  to  which  it  is  to  change.  The  number  appears 
in  5-column  fields  beginning  in  columns  7-11, 
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the  first  four  (7-10)  columns  contain  the 
external  input  number  and  the  fifth  column  (11) 
contains  the  new  value  of  the  output.  This  can 
be  repeated  fifteen  times  per  card;  however, 
only  events  for  the  time  indicated  in  the  first 
five  columns  can  be  on  a card. 

Following  the  circuit  data  cards,  a 7/8/9  card  is 
used  to  separate  the  output  requirements  card,  card  number 
10. 

Card  10:  Element  Outputs  to  be  Plotted.  The  numbers  of 
the  elements  to  be  plotted  appear  in  5-column  fields 
beginning  in  columns  6-10.  Twelve  numbers  are 
allowed  per  card. 

The  final  card  is  a system  control  card.  It  is  a 
6/7/8/9  card. 


Error  Messages 

The  error  messages  can  occur  as  either  fatal  or  non- 
fatal.  Fatal  errors  are  indicated  by  numbers  greater  than 
199.  Non-fatal  errors  are  less  than  199. 

Fatal  Errors 

Fatal  errors  cause  execution  to  stop  in  all  cases 
except  error  203  which  has  an  override  condition.  Normally, 
the  fatal  errors  are  programming  errors.  The  present  fatal 
error  messages  are  as  follows: 


I 

1 

92 

Error  200:  The  user  has  asked  for  an  element  model  that 
I does  not  exist. 

I 

; Error  201:  The  user  has  specified  an  element  that  has  a 

fanout  greater  than  10. 

Error  202:  A memory  type  element  is  out  of  sequence  as 
required  by  data  card  3. 

Error  203:  EVESIM  was  unable  to  initialize  the  circuit. 
Increase  LIMIT  and/or  set  LOVER  equal  to  1 on  data 
card  8. 

Error  204:  Array  LIST (I)  is  full  and  an  attempt  was 
made  to  schedule  another  event.  The  user  has 
reached  the  limit  of  EVESIM.  The  array  LIST (I)  must 
be  increased  in  size. 

Error  205:  The  user  requested  the  output  of  an  element 
that  could  not  be  found  in  the  output  file. 

Non-Fatal  Errors 

Non-fatal  errors  are  errors  that  occur  during  the 
course  of  a simulation  of  a circuit  caused  by  generation  of 
illegal  conditions  within  cr.e  circuit.  These  error  messages 
are  intended  to  aid  the  user  in  identifying  flaws  in  the 
circuit  design  or  implementation  of  a design.  The  existing 
error  messages  are  as  follows: 


Error  1:  The  input  pulse  to  the  clock  input  of  a flip- 
flop  is  too  narrow. 


93 


n 


i 

F 


1 


1 

I 

i 


Error  2:  The  J and/or  K inputs  to  the  JK  master-slave 
flip-flop  are  not  stable  while  the  clock  is  high. 

Error  3:  The  preset  and  clear  inputs  to  a JK  master- 
slave  flip-flop  are  changing  to  a high  value  at  the 
s ame  t ime . 

Error  4:  The  input  pulse  to  the  clear  input  of  a flip- 
flop  is  too  narrow. 

Error  5:  The  input  pulse  to  the  preset  input  of  a flip- 
flop  is  too  narrow. 

Error  6:  The  preset  and  clear  inputs  to  a flip-flop  are 
changing  at  the  same  time  before  either  a preset  or 
clear  input  was  able  to  cause  a change  on  the  output 
of  the  flip-flop. 

Error  7:  The  J and  K inputs  to  a JK  flip-flop  have  the 
same  value.  The  present  state  of  the  flip-flop  is 
unknown;  therefore,  a new  output  cannot  be  deter- 

J 

mined. 


Debug 

During  the  course  of  the  development  of  EVESIM,  a 
debug  ability  was  integrated  into  the  program.  Originally, 
the  intent  for  its  use  was  only  for  development  purposes  and 
was  later  to  be  removed  from  the  program.  However,  it  has 
been  left  in  the  prograr.  af  en  aid  to  continued  maintenance 
of  the  proqram. 


1 
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Debug  basically  consists  of  three  options.  The 
% first  is  to  leave  column  30  blank  on  data  card  1.  This 

effectively  provides  no  debug  information.  A 1 in  column  30 
will  execute  all  portions  of  debug.  If  DEBUG  is  set  equal 
to  2,  array  IDES0U(I,J)  is  not  printed.  All  other  debug 
information  is  printed.  If  DEBUG  equals  3,  array  LCOUNT(I) 
is  printed. 

Figure  26  is  an  example  output  for  DEBUG  equal  to  1. 
The  first  four  lines  is  the  array  IDESCR(I,J)  as  defined  in 

! 

Chapter  3.  The  next  10  lines  is  the  array  IDESOU(I,J)  also 
defined  in  Chapter  3.  The  clock  fanout  is  listed  next.  The 
remaining  information  is  self-explanatory. 

Array  LCOUNT(I)  is  used  as  a bookkeeper.  It  is  used 
to  count  the  number  of  times  a particular  subroutine  is 
called  during  a simulation  execution.  Table  6 defines 
LCOUNT ( I ) . 

Adding  Element  Models 

Table  6 indicates  that  there  are  more  element  models 
available  for  use  than  were  previously  mentioned.  These 
include  a D flip-flop  (DETl) , a 4-bit  counter  (COUN) , a 4- 
bit  parallel-in-parallel-out  shift  register  (SHRE) , a 4-bit 
adder  (ADD4) , and  a BCD-to-decimal  decoder  (DECO) . In 
addition,  an  unused  name  for  a combinational  logic  gate  has 
been  provided  called  UNS.  At  this  time,  these  models  have 
not  been  included  in  EVESIM.  However,  the  overhead  required 
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Table  6.  Definition  of  LCOUNT(I). 


I Subroutine 


1 

Output 

2 

FUNSIM 

3 

JKMl 

4 

DETl 

5 

COUN 

6 

SHRE 

7 

MODI 

8 

SCHED 

9 

UNSCHE 

10 

COMLOG 

11 

NAN  2 

12 

NAN  3 

13 

NAN  4 

14 

NAN  8 

15 

NOR2 

16 

INVE 

17 

XOR2 

18 

OR2 

19 

AND 

20 

UNS 

21 

ADD  4 

22 

DECO 

23 

lODEST 

24 

ERROR 

25 

READ IN 

26 

INCOND 

27 

INITIAL 

28 

PLOTl 
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lias  been  programmed  into  EVESIM  to  include  these  models. 

This  has  been  done  in  order  to  facilitate  the  addition  of 
these  models. 

The  user  of  the  program  is  reminded  that  the  model 
is  required  to  calculate  when  an  event  shall  occur.  The 
actual  scheduling  of  the  event  is  accomplished  through  the 
use  of  subroutine  SCHED.  Certain  information  must  be 
available  to  the  model.  This  information  is  passed  via 
common  Blocks.  The  required  common  Blocks  are  listed  in 
Table  7. 

The  values  of  I0D1-I0D12  have  been  established  by 

subroutine  lODEST.  lODl  is  the  pointer  to  IDESOU(I,J)  of 

the  element  being  processed.  7.002-10012  are  the  pointers  to 

I0ES0U(I,J),  which  are  the  input  elements  to  1001,  listed  in 

order.  The  value  of  the  input  is  found  by  IDES0U(I0D^,  JOS). 

The  new  outputs  are  calculated  and  the  times  of  change  are 

calculated.  The  times  of  change  are  placed  in  lOESOU ( 1001^ , 

J06)  and  the  new  values  are  placed  in  I0ES0UT(I001  , J07) 

n 

and  NOUT(N).  Subroutine  SCHED  is  called  to  schedule  the 
events  in  the  future. 

In  order  to  add  more  models  than  are  listed  in  Table 
6,  additional  changes  must  be  made.  The  size  of  LTYAL(16) 
and  LTN0U(16)  must  be  increased  to  reflect  the  total  number 
of  models.  In  conjunction  with  this  Block  Data  subroutine 
DES  must  be  changed.  LTYAL(n)  must  contain  the  new  element 
name  and  LTNOU(n)  must  ccr.t-.in  the  number  of  outputs  of  the 


( 


Table  7.  Common  blocks  required  by  subroutine  models. 


COMMON/ELEM/LECH,  lETIME 

COMMON/INCL/ICL(100,2) , LIST(1024),  LSP,  LENCL,  INUV(15,2) 
COMMON/MACH/IND , TNI 

COMMON/DESCRP/IDESCR(100,19) , IDESOU (350 , 17 ) , JIl,  JI2,  JI3, 
JI4,  JI5,  JI6,  JI7,  JI8,  J19,  JIIO,  Jill,  JI12,  JI13,  JI14, 
JI15,  JI16,  JI17,  JI18,  JI19,  JOl,  J02,  J03,  J04,  JOS,  J06, 
J07,  JOS,  J09,  JOlO,  JOll,  J012,  J013,  J014,  J015,  J016, 
J017,  LISTCL{50) 

COMMON/WORK/IODl,  IOD2,  IOD3 , IOD4,  IOD5,  IOD6,  IOD7 , IOD8, 
IOD9,  lODlO,  lODll,  IOD12,  II,  NOUT(IO),  NGATE,  lEX,  LFF , 
LOM,  LTYAL{16),  LTNOU(16),  NAVL,  lODL,  NOW,  NSCTI , LERROR, 
NEXIN,  LFATAL 

COMMON /VALE/LOW,  LHIGH,  LUNK,  IFAN,  LIMIT,  LOVER 
COMMON/DEBUG/  LBU,  LCOUNT(40) 

COMMON/WID/NCLWID,  NPLWID,  NCLW2 , NPCW2 


new  element.  NAVL  must  also  be  adjusted  to  reflect  the 
total  number  of  models  available.  With  these  changes, 
additional  models  can  be  added  as  explained  above. 


CHAPTER  8 


[ 
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CONCLUSIONS  AND  RECOMMENDATIONS 

EVESIM  has  been  found  to  be  an  effective  program  for 
analyzing  a digital  logic  circuit.  EVESIM  provides  a means 
for  the  designer  to  quickly  verify  his  logic  circuit  designs. 
EVESIM  also  gives  the  designer  a tool  with  which  his  design 
can  be  optimized  and  quickly  re-verified.  The  time  con- 
suming task  and  manually  tedious  operation  of  design  verifi- 
cation of  logic  circuits  has  been  reduced  to  a computer 
program. 

Conclusions 

The  original  objectives  as  stated  in  Chapter  1 have 
been  accomplished.  EVESIM  is  a level  mode  simulator  that 
uses  element  level  models.  The  efficiency  gained  in  using 
a time-driven  event  simulator  have  allowed  the  original 
objectives  to  be  expanded  to  include  the  use  of  three  valued 
logic  and  the  individual  assignment  of  delays  depicting 
whether  an  element  output  is  rising  or  falling.  In  addi- 
tion, the  use  of  an  iterative  process  during  initialization 
has  allowed  a limited  capability  for  circuit  stability 
analysis . 

The  three  valued  logic  is  used  to  identify  illegal 
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This  information  is  used  by  the  designer  to  improve  his 
design  or  identify  flaws  to  be  corrected  in  his  design.  The 
three  valued  logic  is  also  used  during  the  initialization 
process.  The  initialization  of  the  circuit  can  then  be 
accomplished  within  the  program  without  any  requirements  for 
ordering  of  the  elements  of  the  circuit. 

The  assignment  of  delays  relative  to  whether  a 
signal  is  rising  or  falling  presents  a much  improved  model 
over  other  techniques.  The  modeling  accuracy  gained  is  in 
keeping  with  the  growing  need  for  more  accurate  timing  in- 
formation as  the  speed  of  logic  increases.  This  technique 
also  insures  that  EVESIM  will  not  grow  old  and  useless  as 
loqic  families  change. 

^ Iterations  are  made  through  the  circuit  during 

initialization  until  all  unknown  conditions  are  cleared  and 

no  e'lement  outputs  are  still  toggling,  or  until  the  preset 

leveJ^  of  attempts  to  initialize  is  reached.  If  increasing 

the  npmber  of  iterations  does  not  result  in  the  elimination 
'\ 

of  elements  toggling,  then  those  elements  toggling  are 
flagged  as  being  part  of  an  unstable  condition.  Certainly, 
this  may  be  part  of  the  design,  and  therefore  EVESIM  has  an 
override  capability.  However,  at  a minimum,  the  user  is 
made  aware  of  the  situation  such  that  he  can  make  a decision 
as  to  whether  it  vri.ll  be  e.l?  r'-’ed 

■EVESIM  has  also  pro.'iced  the  user  with  an  improved 


ability  ',to  identify  errors  ’n  >'is  circuit  through  the  use  of 
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error  messages.  In  conjunction  with  this,  debugging  in- 
formation is  available  when  it  is  needed. 

It  is  hoped  that  most  programming  errors  have  been 
eliminated  from  EVESIM.  However,  it  is  a rule  of  program- 
ming that  the  elimination  of  one  error  is  going  to  cause 
another  error  some  other  place  in  the  program.  It  is  time 
that  EVESIM  move  from  a development  phase  to  a phase  of  wide 
usage.  It  is  during  this  phase  of  wide  usage  that  those 
hidden  programming  errors  will  be  identified.  It  is  hoped 
that  EVESIM  will  continue  to  grow  and  be  expanded  to  include 
the  recommendations  of  the  next  section. 

Recommendations 

The  process  of  developing  a program  is  never  ending. 
If  time  were  not  a factor,  the  author  could  continue  forever 
adding  to  and  improving  the  program.  Areas  of  recommended 
improvement  for  EVESIM  are  discussed  here  in  hope  that 
someone  may  be  able  to  include  them  in  a future  version. 

Assignment  of  delay  values  can  be  partially  auto- 
mated. Certainly,  it  would  be  more  realistic  to  assign 
delays  automatically  by  a random  number  generator  process. 
The  limits  of  the  random  number  generator  could  easily  be 
the  minimum  and  maximum  values  of  delays  as  depicted  in  a 
data  book.  This  would  relieve  the  burden  to  the  user  of 
assignment  of  "typical"  delay  values.  In  conjunction  with 
this  process,  the  delay  value  assignment  could  also  be 
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weinhted  by  the  fanout  of  the  element.  This  would  not  be  a 
major  extension  since  EVESIM  already  qenerates  a fanout 
description  for  other  purposes. 

A further  extension  of  the  delay  information  would 

be  to  add  a fourth  logic  va.'.ue.  The  intended  use  of  this 

fourth  value  would  be  to  indicate  that  an  element  output  is 

» 

changing  in  value.  This  could  possibly  be  defined  as  a 
non-propagating  unknown  condition.  It  would  represent  the 
element  output  during  the  time  from  when  an  input  change  has 
occurred  until  the  output  was  stable  at  its  new  value. 

The  timing  mechanism  for  EVESIM  is  no  longer  very 
efficient  if  the  length  of  the  internal  clock  is  forced  to 
become  very  long  as  the  result  of  a single  long  delay.  For 
example,  a circuit  containing  a single  one-shot  element 
would  require  the  clock  to  be  of  the  length  of  the  one-shot. 
This  problem  can  be  overcome  by  establishing  an  events  queue 
for  events  that  are  beyond  the  length  of  the  clocking 
mechanism.  If  these  events  are  placed  in  a stack  and 
ordered  by  time,  then  after  each  pass  through  the  clocking 
mechanism,  the  event  at  the  top  of  the  stack  is  checked  to 
see  whether  the  event  will  fit  into  the  next  pass  through 
the  clocking  mechanism.  If  it  will  fit  in  the  time  frame 
of  the  next  pass,  it  is  then  inserted  appropriately  into  the 
normal  timing  mechanism  for  processing.  This  technique 
v/ould  be  much  more  efficient  than  extending  the  length  of 
the  clock  to  accommodate  a single  element. 
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The  final  recommendations  is  probably  a much  longer 
* range  goal  than  the  previous  recommendations.  The  above 

recommendations  require  minimal  changes  to  EVESIM.  How- 
ever, it  is  considered  that  an  extension  to  include  an 
ability  to  perform  critical  race  analysis  would  be  an 

« 

• t 

! ultimate  and  final  addition  to  EVESIM.  The  inclusion  of 

this  ability  along  with  the  existing  SCRITSS  test  sequence 
generator  and  the  completion  of  the  MACRO  AHPL  compiler 
would  provide  maximum  automation  to  logic  circuit  design. 


) 
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